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1, INTRODUCTION 

Doubly-fed inductor generator (DFIG) wind turbine is the most used wind turbine generator in the 
world today, with its usage rising above 85% of functioning wind turbines [1]. The main advantage of DFIG 
wind turbine is its ability to operate with varying wind speed without stressing its mechanical and electrical 
components. Other advantages of DFIG include; it is less expensive, it 1s capable of controlling reactive power, 
it has high power efficiency and it reduces power variation. Modelling wind energy systems using power 
systems software for the purpose of carrying out several analyses is necessary, considering the high initial cost 
of implementation. Advantages of analyzing power systems on software include and are not limited to the 
following; problems could be predetermined and solutions can be proposed, service time of the equipment can 
be prolonged, response of the system under different scenarios can be evaluated, downtime could be reduced 
and running cost eventually reduced. 

Temporary overvoltage is one of the issues that occur on any power system. It involves oscillation 
which exists for long duration of time, which could be damped or weakly damped. Events leading to temporary 
overvoltage include; resonance, ferroresonance, load rejection, fault clearing and harmonics. Ferroresonance 
is considered in this paper due to the severity of damage on transformers which is one of the most expensive 
equipment in power systems. Ferroresonance could occur on wind power systems due to the interaction of non- 
linear inductor and capacitor when a disturbance occurs in the system. Ferroresonance could lead to high 
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permanent overcurrent and overvoltage, displacement of neutral point voltage, overheating in transformer and 
excessive continuous loud noise from the transformer [2, 3]. Power systems simulation software can be used 
to simulate, analyze potential threats and mitigating solutions can be obtained. 

There are numerous power systems software available today, some of these software could have good 
performance for a particular analysis and may not be effective for other analyses when compared to real power 
systems. ATP/EMTP and MATLAB are packages common among researchers to carry out analyses on power 
systems. Alves and Orizondo [4] compared the use of ATP/EMTP and Matlab/Simulink program to simulate 
unified power flow controller (UPFC). The UPFC was implemented by using two voltage source converters 
coupled through a DC link capacitor. It was found that ATP calculation speed is faster than that of MATLAB 
due to the nature of the solver algorithm. ATP/EMTP and Matlab/Simulink were compared for time domain 
transient studies by Danyek et al. [5]. The two packages were used in lightning surge simulation and it was 
concluded that MATLAB offers more possibilities and options for signal processing and control. Simulation 
of arc faults using ATP and Simulink for comparison was done by Khan and Bengiamin [6]. Cassie arc model 
was simulated in two circuit schematics and the comparison was done in both time and frequency domain. 
Time step of | us and sample time 0.2 s were used for the two packages. The results were similar for both time 
and frequency domain and it was added that the little difference is as a result of the modelling of electrical 
passive components such as capacitors and inductors rather than solver algorithms. Three tools; EMTP, 
MATLAB and XTRANS were compared by Schavemaker, et al. for electrical transient computations. 
Simulation of electric arc using Schwarz arc model with time step of 6 us was performed on all the tools. 
XTRANS was found to be best of the three in terms of user-friendliness, MATLAB was found to be the most 
effective package in terms of programming freedom and flexibility while EMTP gave the best solution and 
could be used for larger simulations [7]. Barsali and Ceraolo [8] did simulation of electromagnetic transient in 
power systems comparing ATP and Simulink using different time steps. It was concluded that simulation took 
less time on ATP while Simulink result is more accurate and simulating large circuit is easier on it. In this 
paper, the two packages are compared for temporary overvoltage, with emphasis on ferroresonance on DFIG 
wind turbine, for accuracy, time of simulation and ease of simulation. Post-analyses were also done for the 
output signals by characterizing the ferroresonance. 


2. ATP/EMTP 

ATP evolved from electromagnetic transient program (EMTP) as experts around the world started 
contributing to it since 1975. ATP was developed for the public in 1984 at Bonneville Power Administration 
(BPA) of Portland, Oregon. ATP has also developed over the years and ATPDraw was birthed at Norwegian 
University of Science and Technology, Norway. ATPDraw is a graphical mouse driven preprocessor, which 
works by dragging and dropping predefined components to design a power network model. ATPDraw converts 
the generated model to input file for ATP simulation [9]. 

ATP simulation is performed as a function of time, using trapezoidal rule of integration to solve the 
differential equations of the power network components. In addition, model verification can be done in 
frequency domain, where users can specify geometry and the specification of the device. Various predefined 
components in ATP include but are not limited to; passive electrical elements such as resistors, inductors and 
capacitors, transformers, rotating machines, switches, electrical sources, surge arresters transmission lines and 
cables. Furthermore, ATP has the capacity to model control systems through the aid of two features; Transient 
Analysis of Control Systems (TACS) and MODELS (a programming language). ATP/EMTP is used for 
various EMTP applications which are not limited to the following; insulation coordination, switching 
overvoltage analysis lightning overvoltage analysis, protective relay modelling, power quality studies and shaft 
torsional oscillation studies. There are post-processors that are interfaced with the ATP which are used to 
display results [9]. 


3. SIMULINK/MATLAB 

Matrix Laboratory which is abbreviated as “Matlab” is technical computing software which started 
approximately 70 years ago based on earlier library such as LINPACK. It started as a program to handle any 
matrix problem but eventually grow up to be used in various mathematical applications. MATLAB as a 
computing environment has other file extensions such as Simulink, Simscape and MUPAD, which makes it 
more robust and enables it to perform more technical computations. Simulink is a graphical environment for 
modelling, simulating and analyzing dynamical systems. Simulink contains preset blocks called libraries which 
contain linear, non-linear, continuous and discrete functions that could be used to represent any system to be 
analyzed. Simulink has gained popularity among researchers working in power systems these days, this is due 
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to its elaborate libraries, ability to customize the libraries and interfacing capability of diverse program 
modules. 


4. RESEARCH METHODOLOGY 

Modelling of DFIG wind turbine and simulation of ferroresonance using ATP and Simulink are the 
methodologies used in this paper. Figure 1 shows the major components for generating power from wind. The 
wind speed must be greater than the cut-in speed (about 12 m/s), so as to create enough torque to run the 
generator. The generator is connected via the gear box which adapts the speed of the turbine to the generator 
speed. The rotor of the generator converted the kinetic energy of the turbine to mechanical energy and the 
electrical energy is obtained from both the stator side of the generator and the rotor side via the back-to-back 
converter. The output power is 0.69 kV which was stepped up to 33 kV by the wind turbine transformer and 


supplied to the grid. 
GEAR 
Ae 


TURBINE 






CONTROL SYSTEM 


Figure 1. DFIG wind turbine configuration. 


4.1. Wind Simulation 

Wind speed at any location varies with time and season. One of the ways to represent wind speed is 
by using probability distribution function (PDF) [10]. The most commonly used PDF is Weibull distribution, 
however, the wind speed can be forecast using Weibull distribution as shown in Equation (1); 
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Where F(V)is frequency of occurrence for wind speed, V is frequency distribution, W is empirical Weibull 
scale factor, K is empirical Weibull shape factor. 


4.2. DFIG Wind Turbine and Generator 
The power output of wind turbine is represented in Equation (2). The DFIG wind turbine generator was 
modelled according to [11]. 


1 
2 


Pn = >pAv3C, (A,B) (2) 


Where A 1s the swept area, p is the density of the air, v 1s the speed of the wind and C, is the power coefficient 
of the turbine, which depends on tip speed ratio, ’ and pitch angle, B. Tip speed ratio can be represented by 
Equation (3): 


) et (3) 
; 2 ; ; 

Where w,. is the angular velocity of the rotor. w, = a and n is revolution/min 

The equation to obtain optimal C, 1s given in Equation (4): 


=f (4) 
CoA B)=Cy (32 — C38 — Cy) + Cod 


Where Ci, C2, C3, C4, Cs and Ce are 0.71, 1.49, 0.56, 13.1, 0.182 and 0.004 respectively. Where, 
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1 1 0.035 (5) 
A, A+0.08B B3+1 





Seta (6) 
P= Ts (7) 
PS 5 ee (8) 


For a lossless generation operating in a steady state; 
Tas (9) 
P.=P.,— P; (10) 
Substituting Equation 7, 8 and 9 into 10, Equation 11 1s obtained. 


(Ce see (11) 
Ws 
Where, P. and P, are power output for rotor and stator of the DFIG respectively, T,, and T,,, are the 
aerodynamic torque of the turbine and electromagnetic torque of the generator, w, is the angular velocity of 
the stator and S is the slip of the generator. 

DFIG is an induction generator having a configuration where its stator is connected directly to the 
load and its rotor, is connected to the load through power converter. The operation of DFIG can easily be 
explained in two axes, d-q coordinate frame. The conversion from ABC to dq is done using park transformation 
in Equation (12). 
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Where Vsa, Vsq, are the stator voltages in d and q axis respectively and Vira, Vra are the rotor voltages in d and q 
axis respectively, R,, R, are the stator and rotor resistance respectively, L,, L, are the stator and rotor inductance 
respectively, L;,, L;, are the stator and rotor self-inductance respectively, L,, is the mutual inductance, Iq, 
I,, are the stator currents in d and q axis respectively I,g, I-g are the rotor currents in d and q axis respectively, 
Psa» Psq are the stator fluxes in d and q axis respectively and @;q, @;g are the rotor fluxes in d and q axis 
respectively. 
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The torque, active power and reactive power generated by the generator can be obtained using 
Equation (23), (24) and (25): 


Lm 
T, = 1.5 Ny is (IsqPsa eee Psq) (23) 


Where, N, is the number of poles of the generator. 


3 
Po, (Veqleget Vea lea) (24) 


3 (25) 
Q.= 5 (Vsqlsa — Vsq ig) 


Table | shows the parameters of the wind turbine DFIG used in both ATP and MATLAB model. 


Table 1. Parameters of DFIG wind turbine 
Parameters Values (Ohm) 
Mutual inductance 0.8436 
Stator resistance 0.0023 
stator inductance 0.029 
Rotor resistance 0.0073 
Rotor inductance 0.033 


The grid parameters used in the simulation 1s shown in Table 2, where the X/R ratio is approximately 
6. However, the parameter of the filter was calculated based on the cut-off frequency of the filter which is 10% 
of the carrier frequency. Hence, the parameter of the filter was calculated based on Equation (26); 


—— (26) 


Table 2. Grid parameters 


Parameters Values 
Grid voltage 33 kV 
Fault level 10 MVA 
Grid resistance 0.2Q 
Grid inductance 4.26 Q 


4.3. Converters and Their Control 

The back-to-back converter 1s rectifier and inverter connected through a DC link. This scheme allows 
control of reactive power at the grid. The rotor side converter (RSC) controls the active and reactive power 
delivered by stator by ensuring that it 1s needed according to the mechanical torque. However, the grid side 
converter (GSC) controls the power delivered to the grid by ensuring the voltage at the DC link is constant. 
The control scheme of both RSC and GSC 1s shown in Figure 2. 
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Where o is the dispersion coefficient of the generator. 
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Figure 2. The control scheme of both RSC and GSC 


4.4. Transformer 

The wind turbine transformer is rated 0.69/33 kV, 2.5 MVA with vector configuration of DY7. 
According to [12], simulation of electrical transients which range in frequencies from dc to some MHz are 
divided into the four categories as shown in Table 3. The analysis in this study could be categorized under low 
frequency transient, thus transformer modelling was done according to the specification of Table 3. 
Nonlinearity characteristic of the transformer is important while modelling for temporary overvoltage, the 
relationship between flux linkage and magnetizing current which can be calculated using Equation (32): 


I = ad+ bA” (32) 


Where / is the magnetizing current, a and b are the constant coefficient for linear and saturated region 
respectively, 4 is the flux linkage and n is number of polynomial degree order of the saturation. 


Table 3. Cigre standard for modelling transformer according to frequency range [12] 


Parameters Low frequency Slow front Transients Fast front transients Very fast transients 
Transients 

Short-circuit impedance Very important Very important Important Negligible 

saturation Very important Very important” Negligible Negligible 

Iron losses Important®) Important Negligible Negligible 

Eddy current Very important Important Negligible Negligible 

Capacitive Coupling Negligible Important Very important Very important 


(1) Only for transformer energization phenomena, otherwise important 
(2) Only for resonance phenomena 


5. MODELLING DFIG WIND TURBINE CONFIGURATION ON ATP/EMTP 

Overview of DFIG wind turbine configuration on ATP is shown in Figure 3, which presents: 
simulation of the wind turbine, DFIG generator, back-to-back controller and control, and wind turbine 
transformer. 
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Figure 3. Overview of ATP model of DFIG wind configuration 


5.1. ATP Model of Wind Turbine and Generator 

DFIG 1s modelled using UM4 tool in ATP. According to ATP rule book [13], when modelling a single 
mass model, the mechanical elements are represented by equivalent electrical variables such that; Moment of 
inertia constant was represented by capacitance, damping coefficient was represented by conductance, while 
aerodynamic and electromechanical torque were represented as current source using (Type 60) and angular 
velocity was represented by voltage. This 1s shown in Figure 4. Where, Ry, 1s small torque sensing resistor, R, 
is large resistor, Rx 1s transmission cable, Dm 1s damping coefficient, Jm is the capacitive element and Ty is a 
current source. 


Current sources 
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& Machine 
UM4 





oe GY 


Figure 4. Electrical representation of mechanical elements of wind turbine and generator 


5.2. ATP Model of Back-to-Back Converter and Control 

Back to back converter which is made up of six pulse rectifier and inverter connected through DC 
capacitor is modelled by the use of TACS-controlled TYPE 13 switch. The operation of the switch is such that 
if one of the switches on a phase leg is opened, the other must be closed at all times. The logic through which 
the switches operate depends on output signal after the comparison of the output voltages of the generator and 
the modulating signal of pulse width modulation (PWM). PWM controls the output voltage for the converters 
by propagating a sawtooth signal which switches off to disallow passage of high frequency voltage output. 
PWM was done using MODEL, where the carrier frequency is 2 kHz. 
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5.3. ATP Model of Transformer 

ATP has three model types for transformer; saturable transformer component (STC) model, BCTRAN 
model and Hybrid model. The BCTRAN model is used to present the transformer. The model can be computed 
by considering steady state single phase multi-winding transformer in a matrix form as shown in Equation (33). 


[V] = [7] [Z] (33) 


The impedance Z in Equation 33, is the both the self-impedance of a phase and the mutual impedance of phases. 
However, for a transient study, Equation (33) would change to (34): 


V] = (1 (R)+ |S] (2) (34) 


Where, R and L are the real and imaginary components of the impedance. 

BCTRAN model can compute the core saturation of a transformer but cannot represent it. Hence, an 
external nonlinear inductance can be connected at the transformer terminal as shown in Figure 5. 

In order to use this model, open and short circuit test data of the transformer is needed. The parameters 
of the wind turbine transformer used for this model is shown in Table 4 as found in [14]. 

The RMS voltage versus current data was converted to resultant A- I data by the BCTRAN supporting 
routine and curve fitting was done to improve the saturation curve as shown in Figure 6. The core saturation 
characteristic is implemented using pseudo-nonlinear hysterestic inductor TYPE 96 of ATP and it was placed 
at the low voltage terminal of the transformer. 


Primary Side Secondary Side 


Equivalent 





Figure 5. Schematic diagram of BCTRAN Model of two windings transformer with external non-linear 
representing the core [9] 
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Figure 6. The nonlinearity characteristic curve of the wind turbine transformer. 


Table 4. The test report of wind turbine transformer DY7 


Test Voltage (%) Current (%) Power loss (kw) 
Open circuit test @ LV winding 90 0.43 aah 
100 0.9 5.5 
110 1.352 8.1 

Impedance (%) Power (MVA) Power loss (kw) 
Short circuit test 8.3 2.5 eal 
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6. RESULTS ON ATP 
6.1. Simulation of Ferroresonance on ATP/EMTP 

Simulation time step of 5 us and sample time 3 s were chosen for the simulation. Figure 7 shows the 
output voltage during normal operation of DFIG wind farm. Similar to [2], two switching events leading to 
ferroresonance were considered; Nuisance operation of the circuit breaker during opening and closing. 
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Figure 7. Output voltage of DFIG during normal operation on ATP 


6.2. Nuisance Operation of the Breaker During Opening and Closing on ATP 

The circuit breaker connecting the secondary side of the transformer was considered in this simulation. 
Either during the opening or closing operation of the breaker, all the three poles are expected to operate, but in 
a situation where one or two pole(s) get stuck during the process, this could lead to ferroresonance. Figure 8a 
shows resulting ferroresonance due to one stuck pole respectively during opening, and Figure 8b shows 
ferroresonance during closing operation. Overvoltage of 4.22 P.U. and 3.77 P.U. was experienced during 
opening and closing event respectively. In both cases, the overvoltage was higher than the power frequency 
withstand voltage of the equipment according to [15] and could be dangerous for the equipment due to the 
undamped nature of the overvoltage. 
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Figure 8. (a). Output voltage due to a stuck pole during opening of circuit breaker on ATP and (b). Output 
voltage due to a stuck pole during closing of circuit breaker on ATP 


6.3. Characterization of Ferrosonance on ATP/EMTP 

Fast Fourier Transform (FFT) analysis was performed to characterize the spectrum of the harmonics. 
Characterization is needed in order to determine the type of ferroresonance existing in the system. FFT analysis 
exists as post-processor tool on ATP, thus FFT was done for the waveform between 0.8 s and 1.5 s and results 
are shown in Figure 9a and 9b for the opening and closing ferroresonant event. In both cases, 1‘, 18 and 20" 
harmonic frequencies are very prominent in the spectrum. Chaotic or quasi-periodic mode ferroresonance may 
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exist in the system due to non-periodic frequency that is continuous. Hence, further analysis is needed to better 
characterize the ferroresonance. 
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Figure 9. (a) FFT analysis of ferroresonant waveform during opening on ATP and (b) FFT analysis of 
ferroresonant waveform during closing on ATP 


7. Modelling DFIG Wind Turbine on Simulink/Matlab 
The model of DFIG wind turbine configuration implemented on simulink is shown in Figure 10. Most 


of the components exist as preset blocks in the simulink library. It entails modelling of DFIG wind turbine and 
generator, converters and their control and three phase, three winding transformers. 









[Vabc_grid_conv] 
[labc_grid_conv] 










Drive Train 








[AsyncMac_sig] 


[Vabc_rotor] 






Pulses_rotor_conv] jPulses__ond_pany} ¢ 





—— 


Breaker 








33kV 





B_rotor_conv Universal Bridge Universal Bridge’ B_grid_conv 





Two winding 
Three phase 
Transformer 






Asynchronous Machine 
pu Units 














Wind Turbine Control 





Figure 10. Overview of DFIG wind turbine configuration on Simulink/MATLAB 
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7.1. DFIG Wind Turbine and Generator Modelling on Simulink 

This was implemented as a block using Equation 1. The wind speed was input by implementing 
Equation 1 and maximum C, is determined by Equation 2 modelled by a look-up wind turbine characteristic 
curve C, (A, B). The output torque was input to the shaft of the generator. The mechanical elements of the 
turbine were also added such that the stiffness of the drive train is infinite, the friction and inertia of the turbine 
were modelled using single-mass model. 

Generator was modelled by using asynchronous machine block already present in the simulink library. 
This could act as either motor or generator depending on polarity of mechanical torque, Tm. Tm was input in to 
this block, while resistance and inductance for both rotor and stator were also input as parameters presented in 


Table 1. 


7.2. Matlab Modelling of Back-to-Back Converter and Control 
Both converters were modelled using universal bridge block and IGBT already exists as a block in 


Simulink. Implementation of converter control as shown in Figure 10 was implemented using Math operation 
blocks. 


7.3. Modelling Wind Turbine Transformer on Simulink 
Wind turbine transformer was modelled using preset three phase, two winding transformer blocks. 
The parameters for the core, primary and secondary winding is shown in Table 5. A nonlinear characteristic 


curve presented in Figure 6 was used for core saturation. 


Table 5. Parameter of wind turbine transformer 


Parameters Primary Secondary 
Vector group D Y 
R (Q) 0.002 0.42 


L (Q) 0.0327 13.41 


Time step and sample time of 5 us and 3 s respectively were chosen for the simulation. The voltage 
during normal operation of DFIG wind turbine on Simulink is shown in Figure 11. 
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Figure 11. Output voltage during normal operation on Simulink 
8. RESULTS ON SIMULINK 
8.1. Simulation of Ferroresonance on Simulink/Matlab 


Figure 12a and 12b show the ferroresonant overvoltage during opening and closing of circuit breaker 
in an event of one stuck pole, with overvoltage of 4.63 P.U. and 6.36 P.U. accordingly. 
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Figure 12. (a) Output voltage due to a stuck pole during circuit breaker opening on Simulink and (b) Output 

voltage due to a stuck pole during circuit breaker closing on Simulink 


8.2. Characterization of Resultant Ferroresonance on Matlab 

Simulink/Matlab can carry out further analysis on the ferroresonance, for the purpose of 
understanding the characteristic of the ferroresonance and to determine the risk level it poses to power 
equipment. Figure 13a and 13b show the FFT analysis performed between 0.8 s to 1.5 s of waveform for the 
opening and closing ferroresonant events respectively. The spectrum of frequency during the opening and 
closing events show that the magnitude are high on 1* , 3% , 5" and 7" order of harmonics. 
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Figure 13. (a) FFT analysis of ferroresonant waveform during opening on Simulink and (b) FFT analysis of 
ferroresonant waveform during closing on Simulink 





However, it 1s quite difficult to characterize ferroresonance by just observing the waveform or by 
performing FFT analysis. Further analyses that could be done are phase plane and Poincaré mapping. Phase 
plane shows the periodicity of signal produced by a non-linear system, since the periodic signal follows a 
closed-loop trajectory [16]. Observation of the trajectory of the phase plane diagram indicates the classification 
of the ferroresonance existing in the system where a single closed path represents existence of fundamental 
ferroresonance mode, while existence of 2 or 3 closed paths represent subharmonic ferroresonance mode. 
Multiple closed paths depict quasi-periodic mode and irregular multiple paths show chaotic mode exists in the 
system. 

Poincaré map is named after Henri Poincaré a French scientist that founded it. It 1s the point of 
intersection of a periodic orbit in signal of a non-linear dynamic theory where a single point of intersection 
depicts fundamental ferroresonance mode while existence of 2 or 3 points represents subharmonic 
ferroresonance mode. Multiple number of points with closed path and multiple points with irregular path depict 
quasi-periodic mode and chaotic mode respectively. Poincaré mapping 1s done by sampling the peaks of voltage 
and output signal at power frequency [17]. Both phase plane and Poincaré mapping were done to further 
characterize the ferroresonance to differentiate between the two modes. Figure 14a and 14b show the phase 
plane diagram for the opening and closing events respectively, while Figure 15a and 15b show the Poincaré 


Indonesian J Elec Eng & Comp Sci, Vol. 14, No. 3, June 2019: 1581 — 1594 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1593 


mapping for opening and closing events accordingly. The ferroresonance mode was chaotic mode for opening 
event because the phase plane diagram had chaotic periodic circular path, which connotes extreme disturbance 
and quasi-periodic mode for closing event because of its closed multiple loops. While the Poincaré map depict 
chaotic mode for opening event since the points are scattered and quasi- periodic for closing event which shows 
multiple number of point clustered together. 
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Figure 14. (a) Phase plane diagram of ferroresonant overvoltage during opening operation on Simulink and 
(b) Phase plane diagram of ferroresonant overvoltage during closing operation on Simulink 
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Figure 15. (a) Poincaré mapping of ferroresonant overvoltage during opening operation on Simulink and (b) 
Poincaré mapping of ferroresonant overvoltage during closing operation on Simulink 


9. CONCLUSION 

There were discrepancies between the two software packages for the two simulated ferroresonance 
events; overvoltage results obtained on ATP were higher than that of Simulink. The time of simulation for ATP 
was 110.58 secs and that of Simulink was 130.81 secs which showed that ATP faster than Simulink. 
Overvoltage of 4.22 P.U. and 3.77 P.U. was experienced during opening and closing event respectively on 
ATP/EMTP model while overvoltage of 6.36 P.U. and 4.63 P.U. was experienced during opening and closing 
event respectively on Simulink model. The difference in results observed between the two packages could be 
from the preset model of the components in the two software and not from the solver algorithm of the software 
this is in agreement with finding of [6]. Any change in parameters when analyzing ferroresonant condition can 
cause change in the behaviour of the system because of extreme sentivity of ferroresonance phenomenon. 
Meanwhile, Simulink was found to be more user-friendly and easier to use for wind turbine simulation due to 
availability of electrical power components as blocks in its library. 

Characterization of the ferroresonance using ATP/EMTP could not easily be determined as FFT 
analysis could only be performed on it. Phase plane diagram and Poincaré mapping performed on Matlab 
showed that chaotic mode and quasi-periodic mode ferroresonance occurred in opening and closing switching 
events rspectively. Thus, Simulink/Matlab allows post-processing of signals for better understanding and 
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characterization of such signals while ATP is limited in this aspect. However, signal output from ATP can be 
imported to Matlab for further analysis. 
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